LA-UR-94-2124 

Chaos and noise in a truncated Toda potential 



Salman Habib*, Henry E. Kandrup^^+, and M. Elaine Mahon 

*T-6, Theoretical Astrophysics 
and 

T-8, Elementary Particles and Field Theory 
Los Alamos National Laboratory 
Los Alamos, New Mexico 87545 

^Department of Astronomy and Institute for Fundamental Theory 

University of Florida 
Gainesville, Florida 32611 



Abstract 

Results are reported from a numerical investigation of orbits in a truncated Toda 
potential which is perturbed by weak friction and noise. Two significant conclusions 
are shown to emerge: (1) Despite other nontrivial behaviour, configuration, velocity, 
and energy space moments associated with these perturbations exhibit a simple scal- 
ing in the amplitude of the friction and noise. (2) Even very weak friction and noise 
can induce an extrinsic diffusion through cantori on a time scale much shorter than 
that associated with intrinsic diffusion in the unperturbed system. 
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In the past several years, substantial interest has centered on the study of stochas- 
tic orbits in nonintegrable Hamiltonian systems. In particular, much work has focused 
on the phenomenon of intrinsic diffusion through cantori, e.g., using the "turnstile" 
model of MacKay et al [I], which leads to interesting scaling behaviour [2], and, more 
recently [3], through the utilisation of local Lyapunov exponents [4]. However, with 
a few notable exceptions (cf. [5]), little has been done to determine how the results 
derived from such analyses are changed if non-Hamiltonian perturbations are allowed. 
This is an issue of key importance given the modern perspective on the connection 
between dynamical chaos and the foundations of classical statistical mechanics [6]. 

In modeling a physical system in terms of a few degree of freedom Hamiltonian, one 
is usually implementing a reduced description, e.g., in terms of collective coordinates, 
that neglects (hopefully) weak couplings to an external environment and/or various 
(relatively) unimportant degrees of freedom. However, such weak corrections can 
oftentimes be viewed as a source of friction and noise, related via a fluctuation- 
dissipation theorem, and, as such, it is important to investigate the qualitative and 
quantitative effects that arise if the Hamiltonian model be perturbed to allow for weak 
friction and noise. Naively, one might expect that weak coupling to an environment 
should only have effects on very long time scales. However, as will be shown here, 
for nonintegrable systems such effects can in fact be important on time scales much 
shorter than the natural relaxation time scale. 

The system considered here in detail is the sixth order truncation of the Toda [7] 
potential, which leads to the Hamiltonian 

H = \[vl + vf) + + y^) + x^y - + + x^y^ + \y^ 

+ X^y + - ^ 1^6 ^ ^4^2 ^ 1 ^2^4 ^ 11^6^ 

where {x^v^} and {y,Vy} represent conjugate pairs. The Hamiltonian evolution was 
perturbed by allowing for a constant friction — ?yv and delta-correlated additive white 
noise F(t), related via a fluctuation-dissipation theorem with temperature ~ 
Other forms of friction and noise may give different conclusions, a possibility currently 
under investigation. 

For each initial condition, an unperturbed Hamiltonian trajectory was computed. 
Multiple Langevin simulations were then performed, using a numerical algorithm (cf. 
[8]) which generates a random F with the proper first and second moments. A time 
step h = 0.001 was used for most simulations and it was verified that smaller time 
steps led to no statistically significant differences. The total time for each integration 
was t < 140 (the dynamical or crossing time is of order unity). 
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Viewed over long time scales, the unperturbed Hamiltonian trajectories divide 
naturally into only two classes, namely regular orbits, with vanishing Lyapunov ex- 
ponent, and stochastic orbits, with nonvanishing Lyapunov exponent. However, on 
shorter time scales, the stochastic orbits divide in turn into two relatively distinct 
types, namely filling stochastic orbits which travel unimpeded throughout the stochas- 
tic regions and confined, or sticky, stochastic orbits, which are trapped near islands 
of regularity by cantori, and only escape over much longer time scales. It is there- 
fore meaningful to consider the effects of friction and noise separately on these three 
different orbit classes. 

Experiments focused on the energy range 10 < < 100, with I0~^^ < ?7 < 10""^ 
and O.I < Q I E < 10. 0. Two different types of experiments were performed: 
(I) For each of a large number of individual initial conditions, corresponding to both 
regular and stochastic orbits, = 48 different noisy realisations were effected. The 
orbits were compared to the unperturbed trajectories and analysed statistically to 
extract the first and second moments of such quantities as the position, velocity, and 
energy, e.g., 



(2) The near-invariant distribution associated with the Hamiltonian evolution was 
sampled to extract 400 initial conditions corresponding to stochastic orbits, and 100 
noisy realisations were effected for each of these initial conditions. The outputs were 
then analysed in two ways, (a) by comparing with the unperturbed trajectories to 
extract first and second moments which average over the ensemble of initial conditions, 
and (b) by binning the orbital data at fixed time intervals to study systematic changes 
in the form of the near-invariant distribution. This latter set of experiments was 
performed for = 30 and 75 with Q = E and rj = I0~^, I0~^, and 10""*. A single 
experiment was done with £" = = 10 and rj = I0~^. Each of these latter experiments 
took approximately 250 node hours on the Los Alamos CM5. 

Viewed in energy space, weak friction and noise serve to induce a classical diffusion 
process. Specifically, at least for early times, when SErms/ Eunp ^ I, SE^ms satisfies 
the simple scaling relation 



where A[E) is only weakly dependent on E. This scaling holds for all three classes 
of orbits, regular, sticky stochastic, and filling stochastic. Moreover, it holds for 
individual initial conditions as well as for ensembles of initial conditions. Viewed in 
energy space, one cannot distinguish between different orbit classes. 

While for individual initial conditions, A exhibits some variability (~ 10%), when 





SE^l^ = A'{E)E^Qt, 



(3) 
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averaged over the ensembles of initial conditions, the best fit value is very well de- 
termined. For E = 75, 1.68 < A < 1.71, for E = 30, 1.66 < A < 1.69, and, for 
E = 10, 1.66 < A < 1.67. There is no systematic residual dependence on rj. The 
results suggest that, when E is decreased, the best fit value of A also decreases, al- 
beit very slowly. In the limit — )• 0, orbits essentially move in a harmonic oscillator 
potential, and one would anticipate an asymptotic approach to the exact oscillator 
result, SE^^^ = 2Ei]Qt^ i.e., A = \/2- (This result holds in the limit ?y ^ 1, assuming 
"orbit averaging," so that, e.^., sin^ t — )• 1/2.) The data are consistent with this 
expectation. 

Viewed in configuration or velocity space, friction and noise have more complicated 
effects, the form of which depend on the class of orbit, i.e., whether regular, sticky 
stochastic, or filling stochastic. For regular orbits the second moments in position and 
velocity grow as a power law in t, albeit more rapidly than the analytically predicted 
relation Sx^ms^ Sv^^rms^ •••oc i^^^ which is satisfied by the integrable cases of a harmonic 
oscillator or a free particle. By contrast, for stochastic orbits these rms quantities 
grow exponentially at a rate A that is comparable in magnitude to the Lyapunov 
exponent x, which characterises the average instability of the Hamiltonian orbit (See 
Fig. 1). This conclusion holds both for entire ensembles and for individual initial 
conditions. In the latter case, one also observes a direct correlation between the 
growth rate for an individual initial condition and the local Lyapunov exponent Xi 
for that initial condition. This correlation is particularly strong for the case of weak 
friction and noise. Moments for the confined stochastic orbits exhibit an intermediate 
behavior. 

Despite these differences, one observes a universal characteristic relating the mo- 
ments for all three classes of orbits, namely a simple scaling in terms of and rj. 
Specifically, provided that the deviations have not become macroscopic [i.e., assum- 
ing SxrmsT Syrms ^ 1), all three types of orbits satisfy a scaling relation 

^X^Yns^ ^Vrms^ ^^x,rms^ ^^y,rms OC E[E,t^, (4) 

where a = b = 0.50ib0.01 for the range of values probed. This was confirmed both 
for multiple realisations of individual initial conditions and for ensembles of initial 
conditions corresponding to stochastic orbits. It follows that, even for stochastic 
orbits, one observes the same scaling in and rj as for regular orbits in a harmonic 
oscillatory potential, although the time dependence is extremely different. 

Another common feature of all three orbit types is that, at sufficiently early times, 
the time dependence is reasonably well fit by a power law, i.e., 

^X^Yns^ ^Vrms^ ^^x,rms^ ^^y,rms ^ E[E^Q Tj t . (5) 
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The interval over which such a fit is appropriate depends on the orbit type: fiUing 
stochastic orbits rapidly begin to show an exponential divergence, confined stochastic 
orbits only somewhat later, and regular orbits continue to manifest an approximate 
power law growth. The best fit value of c depends on the sampling interval: for all 
types of orbits, if one fits to longer time intervals (again before the perturbation has 
become macroscopic), the best fit value of c increases, even for the case of regular 
orbits. Fitting for relatively short times yields a value c 1.10 — 1.25. Fitting regular 
orbits over longer times yields c 1.20 — 1.50. 

To test whether the early time scaling might be a numerical artifact, e.g., due to 
the finite step size, simulations were performed for a harmonic oscillator potential, 
selecting = = 30 and rj = 10~^. The second moments in position and velocity 
were found to be in complete agreement with the analytic predictions for the moments 

^'^rms^ ^Vrms^ 8^x,rms^ and ^Vy^rms- 

To investigate the generality of these scaling relations, simulations were also per- 
formed for the dihedral DA potential of Ref. [9], with a choice of parameter values 
which leads to a large amount of stochasticity. Specifically, the Hamiltonian was 
taken to be of the form 

TT -'- / 2 I 2\ 2 I 2\ ^/ 2 I 2\2 ^22 

H = :^{v, + Vy) - y(x +y ) - j{x +y ) - -X y , (6) 

with M = 2, a = -1, and b = 1/2. 

Analysis of these simulations shows once again a scaling of the form given by (3) 
and (4), with the same values a = b = 0.50ib0.01. The best fit c for Fqn. (4) shows 
more variability from run to run, so that a detailed comparison is difficult to effect. 
However, the data are consistent with the conclusion that c is the same for both 
models. 

When viewed in energy space, friction and noise again induce a simple diffusion 
process of the form (1) for the DA potential. However, in this case the value of A[E) 
shows somewhat more variability. For 1 < < 50, A varies from ~ 1.6 — 2.6, with 
the largest values correlating with energies where, as characterised by the magnitude 
of the Lyapunov exponent, the phase space is most chaotic. 

Generic ensembles of initial conditions corresponding to filling stochastic orbits, 
when evolved with the truncated Toda Hamiltonian (1), exhibit a coarse-grained 
exponential evolution towards an approximately invariant distribution [3]. For high 
energies E > 50, where the regular regions are very small, this distribution appears to 
correspond to a true invariant measure. However, for lower energies, where the relative 
measure of the regular regions is larger, this distribution slowly changes in form over 
longer time scales. As a result of intrinsic diffusion, the orbits can pass through the 
cantori to occupy phase space regions near the regular islands from which they were 
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excluded at earlier times. However, friction and noise can significantly accelerate 
the changes in this near-invariant distribution by serving as a source of extrinsic 
diffusion. This has important implications for problems in galactic dynamics, which 
will be discussed elsewhere [fO]. 

The Hamiltonian evolution conserves energy, restricting orbits to a constant en- 
ergy hypersurface. When friction and noise are included, the energy is no longer 
conserved. However, it is approximately conserved over sufficiently short time scales, 
so that one can still speak of an "almost constant energy hypersurface." It is therefore 
meaningful to quantify the degree to which friction and noise alter the form of the 
near-invariant distribution on a near-constant energy hypersurface. 

At a coarse-grained level, an ensemble of orbits can be characterised numerically by 
binned, projected distributions, such as /(x, j/, t; At) or /(j/, Vy^ t; At), constructed by 
averaging the binned orbital data over some time interval At [3]. In order to compare 
two different distributions fi and /2, one requires a notion of distance. This was 
provided through the introduction of a coarse-grained norm: For two identically 
normalised distributions, y) and f2{x^ j/), binned in an n X n grid of cells of size 



It is with respect to this measure of distance that the Hamiltonian flow evidences an 
evolution towards a near-invariant measure. 

However, if this near-invariant measure be evolved into the future, allowing for 
even weak friction and noise, one can observe systematic changes on time scales much 
shorter than the intrinsic diffusion timescale. The friction and noise can induce a 
significant extrinsic diffusion which allows filling stochastic orbits to become confined, 
thereby populating regions of the phase space near regular islands which were avoided 
by the deterministic near-invariant measure. 

For the Toda potential, the observed changes are larger for lower energies, where 
the regular regions comprise a bigger fraction of the available phase space, and they 
occur more rapidly when the friction and noise are larger in amplitude. The case of 
E = 30, where cantori are very important, is particularly illuminating. For rj = I0~^, 
there is clear evidence that the noisy ensemble is evolving towards a new distribution 
which, over time scales t ~ 100, is approximately time-independent. For rj = I0~^, 
the ensemble again evidences a non-trivial time evolution, but the changes are suf- 
ficiently slow that one does not see clear evidence of an approach towards a new 
time-independent form. For i] = 10""*, changes in the energy are so large that it no 
longer makes sense to speak of an approximately constant energy hypersurface. How- 
ever, the data are still consistent with an evolution towards a modified distribution 



{Ax, Ay}, 
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in which cantori around the regular regions have been breached. 

This behaviour is illustrated for = 30 in Fig. 2. The first panel shows the form 
of the near-invariant distribution fo{x^ y) associated with the deterministic evolution. 
The second panel shows the form of the noisy near-invariant distribution fn{x, y) for 
J] = 10~^. The final panel exhibits the difference fo[x^y) — The determin- 

istic distribution fo[x^y) has four sharp relative minima which have become blurred 
somewhat in the noisy fr^{x^ y) by the diffusion of orbits into lower density regions. 

At least for E = 30, the evolution from /o to is well fit by an exponential, with 
a rate A that scales roughly as log rj. For a f X f binning, the best fit values for the 
slope are: for r] = 10'^, A = -0.0f43; for i] = 10'^, A = -0.0339; and for i] = 10'^, 
A = -0.0460. 

The fact that orbits can pass through cantori, going either in or out, implies the 
possibility of changes in orbit class between filling and confined stochastic orbits. It 
is difficult to construct a simple numerical algorithm to decide when an orbit has 
changed class. However, visual inspection of a large number of orbits (~ lO"*) leads 
to qualitative conclusions that are easily summarised. At high energies, E > 50, 
changes in orbit class are infrequent since the size of the phase space region restricted 
by cantori is relatively small. However, for lower energies, these regions become larger 
and changes in orbit class more common. 

This can be quantified by determining the minimum amplitude of friction and 
noise required to induce one or more changes in orbit class for a significant fraction 
of the orbits within time t = 100. Consider, e.^., the case E = 20 and ~ For 
T] much less than I0~^, the friction and noise are too weak to cause a significant 
number of orbits to change class. However, for rj ~ I0~^, friction and noise begin to 
become more important, and, already for rj as large as I0~^, as many as 50% of the 
noisy realisations for any initial condition can result in a change between filling and 
confined stochastic. Such changes are not accompanied by changes between regular 
and stochastic orbits, which, deterministically, are separated by KAM tori, rather 
than cantori. Only for rj > 10""^ are any such changes observed. 

What these simulations imply is that even very weak friction and noise, with 
characteristic time scale tji > 10^ characteristic crossing times tcr can have significant 
effects within a time as short as ~ 10^ t™. 
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Figure 1: The early time exponential growth of Sr^ms = {^^Ims + ^vlmsY^"^ (dashed 
curve) and the corresponding Sv^ms (solid curve) for = 30 and rj = 10~^ is clearly 
apparent in this figure. 
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Figure 2: (a) A grey scale plot of the deterministic near-invariant fo{x^ y) for E = 30, 
generated from a 40 X 40 binning. Darker shades represent higher densities. 
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